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I. INTRODUCTION 


Unquenched simulations with fermions that satisfy the Ginsparg-Wilson relation [l| are 


ion [l| 


computationally much more demanding than those with staggered or Wilson fermions. A 
precise comparison between the numerical cost depends on many details of the simulation, for 
example the lattice spacing, the quark masses and also the practical implementation of the 
Gin^arg-Wilson (GW) fermion, he., how well the overlap RES or the domain-wall fermion 


3 S is 


is approximated. A recent review E) however, suggests that dynamical Ginsparg- 
Wilson fermions may be about ten to hundred times more expensive than corresponding 
simulations at comparable masses with either improved staggered fermions using the Asqtad 


action 


10| or twisted-mass Wilson fermions 


0 . 


For that reason most dynamical 


GW simulations so far have been carried out on small volumes together with rather heavy 


quark masses (see Ref. E references therein). Simulating volumes and quark masses 


comparable to those in present-day staggered simulations [l^, for example, is out of reach 
in the near future. 

A computationally affordable compromise for certain applications may be so-called mixed 
fermion simulations. This type of simulation employs GW fermions only for the valence 
quarks, while the sea quarks are either staggered or Wilson fermions. In such an approach 
at least the valence sector exhibits all the benehts stemming from the exact chiral symmetry 


at non-zero lattice spacing 


12j | . Moreover, provided one can use already existing unquenched 


conhgurations generated with either staggered or Wilson fermions, mixed simulations only 
require the computation of correlators in the background of these conhgurations. The ad¬ 
ditional numerical cost is therefore comparable to quenched GW fermion simulations. Such 


mixed actions in two dimensions were studied in the Schwinger model in Ref. 


In four 


dimensions, preliminary results using the publicly available MILG conhgurations tog ether 




with domain-wall or overlap valence fermions have been reported recently 

In this paper we construct the low-energy chiral ehective Lagrangian for a mixed lattice 
theory with staggered sea and GW valence quarks. Based on this ehective Lagrangian 
we compute the pseudoscalar meson masses and decay constants to one loop. Our results 
provide the leading quark mass and lattice spacing dependence of these quantities. The 
formulae are needed in order to analyze numerical data of mixed simulations. 



This paper parallels Refs. [121, llSl], where the chiral Lagrangian for the mixed lattice 
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theory with Wilson sea quarks and GW valence quarks was constructed. We first find all 


the relevant operators of O(a^) in the Symanzik effective action 
theory. A spurion analysis similar to that in ordinary XPT 


21 


igjMfor 1 

L 22, 3 i 


for the mixed lattice 
is applied to the 


Symanzik effective action j24l.l25|. The result is a chiral effective theory that exhibits explicit 
quark mass and lattice spacing dependence of the underlying lattice theory. A recent review 
of this method can be found in Ref. Q. 

The leading order chiral Lagrangian presented here includes the lattice spacing effects 
proportional to a?. Compared to the leading order Lagrangian of staggered XPT 0 0 



1291. |2i|, which is the appropriate low-energy effective theory for “unmixed” lattice QCD 
with staggered sea and valence quarks, it contains only one additional operator together 
with an undetermined low-energy constant. This new operator contributes to the masses 
of “mixed” mesons (one valence and one sea quark) at tree level in the chiral expansion. 
The decay constant of valence-valence mesons then receives contributions from the new 
operator through the mixed meson masses at one loop. However, the masses of valence- 
valence mesons themselves get no one-loop contribution from the new operator. Besides an 
analytic contribution of the one-loop pseudoscalar masses depend only on the 

leading order low-energy constants of staggered XPT and the next-to-leading order (NLO) 
low-energy constants of continuum XPT, the Gasser-Leutwyler coefficients. To the extent 
that these low energy constants are known from previously performed staggered simulations 
0|, the quark mass and lattice spacing dependence of the valence-valence pseudoscalar 
mesons in the mixed theory are highly restricted and depend on only one free parameter. 

The paper is organized as follows: In section E] we discuss the Symanzik effective action 
for the mixed lattice theory and list all those operators that are relevant in the subsequent 
analysis. We then perform the necessary spurion analysis and derive the O(a^) terms in the 
chiral effective Lagrangian. The tree-level meson masses are written down in section IIII A1 
In section IIIIRI we calculate the valence-valence meson masses at one loop; while the cor¬ 
responding calculation of the pseudoscalar decay constants is given in section IIII Gl We 
conclude with a discussion of the results in section IIVI Appendix is devoted to some 
technical details of our calculation. 
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II. THE CHIRAL EFFECTIVE LAGRANGIAN 


A. Lattice theory 


Mixed fermion theories are a generalization of partially quenched theories. Theoretically 
they can be formulated b y a n action with anti-commuting sea and valence quarks and 
commuting valence ghosts j3l|, where the quark masses and the Dirac operators are chosen 
differently in the sea and in the valence sector.^ In the following we consider a mixed theory 
with Nf staggered sea and Ny valence fermions. 

The sea sector is described by the standard staggered fermion action, 


5'sea = +mses.]xsix), 


( 1 ) 




where denotes the gauge covariant central difference operator and msea is the Nf x 
Nf mass matrix in the sea sector. For brevity we have suppressed the flavor and color 
indices. This action is invariant under lattice rotations, axis reversal, and translations by 
one lattice spacing (the so-called shift symmetry). In addition, the single-flavor staggered 
theory possesses vector and (in the massless case) axial-vector U{1) symmetries. The explicit 
expressions for the held transformations that correspond to these symmetries can be found 
in Ref. 2. For Nf havors, these symmetries extend to a U{Nf)i x U{Nf)^ symmetry 
that corresponds to havor transformations on the odd and even sites separately.^ 

The action for the valence and ghost quarks is given by 


5'val — o!^ ^ '4’vi^) {^GW + "^Val (l “ \o.Dgw) } . 


( 2 ) 


The valence and ghost quark masses are contained in the 2Ny x 2Ny mass matrix ruvai- The 

operator that satishes the Ginsparg-Wilson 


relation Q, realized by either overlap |^, l2l, 1^ or domain-wall fermions [^. Crucial is that 
the valence action has an exact chiral symmetry if the valence mass is set to zero |^. In 
addition it is invariant under the lattice symmetries (translations, rotations and rehections). 


a 

oca 

B 

2,^ 




^ Instead of this “ghost” formulation one could also employ the “replica method” with valence quarks only 

Q 

^ The subscripts £ and r are used instead of L and R because the chiral rotations in U{Nf)i x U{Nf)r act 
on the spin and taste degrees of freedom (see Ref. 
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B. Symanzik action 


At momenta p much below the lattice cut-off momentum 1/a, the long dista nce phy sics of 


the lattice theory can be described by the continuum Symanzik effective theory 


19 


The 


effects due to a non-zero lattice spacing appear in the form of higher dimensional operators in 
the effective action, multiplied by appropriate powers of a. These operators are constrained 
by the symmetries of the underlying lattice theory. The Symanzik effective action for the 
mixed lattice theory has the generic form 


‘S'sym — S4, + a Sq + ■ ■ ■ 


( 3 ) 


The hrst term is the known continuum partially quenched QCD action |31l. l34| containing 
4iVj sea quark helds -^g, -03, and Ny valence quark and ghost helds 'ipy, Each staggered 
flavor held comes in four different tastes, hence the four-fold degeneracy in the sea quark 
sector. The mass matrix consists of the renormalized quark masses proportional to the 
bare lattice quark masses. The symmetries of the lattice action forbid any dimension three 
operator that would lead to an additive mass renormalization. 

There are no terms linear in a in eq. because no dimension hve operator is compatible 
with the symmetries of the underlying lattice theory. Dimension hve quark bilinears with 
two valence helds are ruled out by the chiral symmetry in the valence sector l3.^ . and 
staggered quark bilinears are forbidden by the axial U(l) and the shift symmetry [dfj. Q- 
Mixed bilinears with one sea and one valence held are not compatible with the separate 
havor symmetries in the sea and valence sector. 

In order to hnd the terms in S'g it will be convenient to distinguish three types of operators: 
Operators that involve only sea quark helds, operators that contain only valence quark 
helds and those that contain both.^ The terms of the hrst two types have been constructed 
previously and can be found in the literature. The operators involving only sea quark helds 
are listed in Refs. ^ for the Nf = 1 case, and the results were generalized to the 


arbitrary Nf case in Refs. |39|. Similarly, the operators containing only valence helds 


^ We collect a valence quark field tpy ^>^3 the associated ghost field 0®^ in one valence field = {tpy, 
and similarly for the anti-valence fields. 

^ We do not need to discuss purely gluonic operators, which also appear at O(a^), since they transform triv¬ 
ially under chiral transformations and therefore will not affect the form of the chiral effective Lagrangian. 
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f. Q, 


are listed in Ref. [18||, where the Symanzik action for the mixed lattice theory with Wilson 
sea and GW valence quarks was constructed. 

What remains to be done here is the construction of the mixed operators containing sea 
and valence helds. Bilinears of dimension six with one sea and one valence held are ruled 
out by the separate havor symmetries in the sea and valence sector, analogously to the 
dimension hve bilinears. The mixed operators are therefore four-fermion operators that are 
products of two bilinears, one from the sea and one from the valence sector. 

We construct these four-fermion operators by closely following the procedure described 


f. 0- 


in Ref. |23. We hrst construct all relevant lattice operators that are compatible with the 
symmetries of the lattice theory and correspond to dimension 6 four-fermion operators in 
the continuum limit. Sending then a to zero gives the desired terms in the Symanzik action. 
This procedure was used in Ref. 0] to construct all four-fermion operators made of two 
staggered quark bilinears. The method is easily adapted to the mixed operators. We present 
the details of the construction in Appendix El and quote here the hnal result. 

The allowed mixed four-fermion operators are of the form 

C'mL = ^s(7Spin®tcoior)V'S?v(7Spin®tcoior)^V- (4) 


The matrix 'jspin represents one of the sixteen Clifford algebra elements and denotes a 
color gauge group generator.® The Dirac and color indices are contracted in such a way that 
the four-fermion operator is a singlet under 0(4) rotation symmetry and color. Furthermore, 
in writing eq. 0 it is understood that the bilinear 'il^y{'jspin G) thoiJtv is an SU{Nv\Nv) 
flavor singlet and that the bilinear 'ip^ilSpin ® tcoior)'^s is an SU{4:Nf) taste singlet. It is 
worth emphasizing that the mixed four-fermion operators do not break the taste symmetry 
in the sea quark sector. 

The separate axial symmetries in the sea and valence sector imply that the bilinears in 
eq. 0 transform either as a vector or an axial vector. Writing lSpin®i^oior laore compactly 
as 75 ^“ we therefore hnd only four mixed four-fermion operators that are allowed by the 
symmetries: 

OmL,! = (^s7/.^5)(^y7/.^v'), oSL,3 = (5) 

OmL,2 = (V's7/.75^5)(V^y7/.75V’y), OSL.4 = (V's7/.75i“^s)(V'y7M75^“V’v), 


® The identity in color space, for which we use the notation is also allowed here. 
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where we have explicitly separated out bilinears containing the color identity matrix, and 
now restrict to be traceless, summing over a. 

C. Spurion analysis 

The leading term in the Symanzik action, 5*4, is just the continuum action of partially 
quenched QCD. In the massless limit it is invariant under the flavor symmetry group 

GpQQCD = SU{ANf + Nv\Nv)l ® SU{ANf + Nv\Nv)r, (6) 

which is expected to be spontaneously broken to its vector part SU{4:Nf + Nv\Nv)v- The 
low-energy physics is therefore dominated by Nambu-Goldstone bosons. These pseudoscalar 
bosons acquire small masses due to non-vanishing quark masses and a non-zero lattice spac¬ 
ing. The latter contribution has its origin in chiral symmetry breaking terms in Sq. 

To construct the chiral Lagrangian that describes these pseudoscalar bosons we follow the 
standard procedure of a spurion analysis. The coefficient Cj of each term Oi in the Symanzik 
effective action is promoted to a spurion held that transforms under havor transformations 
in eq. in such a way that the product CjOj is invariant. The chiral effective Lagrangian is 
constructed from the pseudoscalar helds and the spurion helds with the requirement that it is 
invariant under havor rotations. Once the chiral Lagrangian is constructed the spurion held 
is set to its original constant value. This guarantees that the chiral Lagrangian explicitly 
breaks the chiral symmetries in the same manner as the underlying Symanzik ehective action, 
and reproduces the same Ward identities. 

In order to perform the spurion analysis for the mixed theory it is convenient to introduce 
the following notation. We collect the quark and anti quark helds in single helds, 

T = (V’5,^y), T = (7) 

where tl^v contains both the anticommuting valence quarks and the commuting ghost helds. 
The mass matrix is given by m = diag(ms', my), with ms being the 4iVj x ANf mass matrix 
for the sea quarks and my is the 2Ny x 2Ny mass matrix for the valence quarks and valence 
ghosts. We also introduce the projection operators 

Ps = diag(/ 5 ,0), Py = diag(0, ly) , (8) 
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where Is denotes the ANf x ANf identity matrix in the sea sector, and Jy the 2iVy x 2iVy 
identity matrix in the space of valence helds. 

For our purposes it will not be necessary to construct all spurion helds that make the 
Symanzik action in eq. m invariant. Most of the analysis has already been done and we can 
rely on previously published results. The results of the spurion anal ysis to 0{a^) for the case 
with staggered sea and valence quarks were written down in Ref. j29|; the analysis can be 


found in detail in Ref. 


39|, which also works to higher order. Since now only the sea sector 


contains staggered quarks all we need to do is to include the projector Ps appropriately in 
all spurion helds associated with sea quarks in this reference. These spurion helds render 
invariant all terms in the Symanzik ehective action that are built only of sea quarks. To 
illustrate this point consider the mass term for the sea quarks, i/jgms'ips- In Ref. jSSl this term 
is made invariant by promoting the mass to a spurion held that transforms asms ^ LmsR^ 
under left- and right transformations L and R. In order to make use of this spurion held 
in the mixed theory we write 'ips^s'4^s = Ps'^Ps^ and assume the same transformation 
property as before, he., m —>■ LmRK The constant value to which the spurion is assigned in 
the end, however, is now PsmPs- One can proceed analogously with all the other spurion 


helds in Ref. 


39| |. From the results in Ref. |l^ for the mixed theory with Wilson sea quarks 


and GW valence quarks we can directly determine the spurion helds that are necessary to 
make the valence held operators invariant. What we need in addition are the new spurion 
helds that make the mixed four-fermion operators in eq. o invariant. 

We want to emphasize that the mixed four-fermion operators do break the symmetry 
group GpQQCD, even though each bilinear in them couples helds of the same chirality only. 
The reason is that any four-fermion operator that is invariant under all transformations in 
GpQQCD niust be of the form 


{mf = {^Psr^PsY + {^l^vr^l^vY + 2{^Psr^l^s)i^Pvr^l^v), 


(9) 


where F represents one of the four combinations 75 ^“ in eq. ©• All three types of four- 
fermion operators on the right hand side of eq. Q appear in the Symanzik ehective action. 
However, because the lattice theory does not posses any symmetries relating the staggered 
sea and the GW valence fermions, they do not enter with the same coefficient in front in 
order to sum up to the square on the left hand side. Gonsequently, although an arbitrary 
linear combination of the three operator types is invariant under the subgroup SU{ANf)L ® 





SU{4:Nf)ji ® SU{Nv\Nv)l ® SU{Nv\Nv)r of Gpqqcd, it is not invariant under the larger 
symmetry transformations of (jpoo rn t hat mix the sea and valence sector. 

Following the notation in Ref. 18|, the mixed four-fermion operators can be made 
invariant under arbitrary chiral flavor transformations L G SU{4Nf + Nv\Nv)l and 
R G SU{4:Nf + Nv\Nv)r by introducing the spurion helds 

D = Di® D 2 ^ LDiL^ ® LD2L\ E = El® E 2 ^ REiR^^ ® RE2R\ 


E = El 


Dn — En — En — Gn — 0? P > 


LEiL^ ® RE2R\ 

Pv. 


G = Gi®G 2^ RGiR^ ® LG2L\ 


( 10 ) 


The constant values to which the spurion helds are assigned to in the end carry the 
subscript “0”. The spurion helds transform as 4-tensors under chiral havor transforma¬ 
tions and therefore carry four indices, which need to be properly contracted with the 
indices of the fermion helds in order to form invariants.® For example, decomposing 
C’mL.i = (' 057 /^V’ 5 )(^y 7 M^y) = (^7/.'0)s(^7MV')y in chiral components we obtain 

+ (V’ij7/.V’i?)s(^/’L7/.V’L)y + {'ipRl^.'ipR)s{'ipRl^.'ipR)v■ 


Both bilinears in the hrst term on the right hand side couple left-handed helds only. It is 
made invariant with the spurion held D, where the indices are contracted according to 


D{'^pLl^,'^pL)s{'tpL'l^^i’L)v = {'tpLDi-^f,'4)L)s{'tp LD2'^^^,'^pL)v ■ ( 12 ) 

The other three terms in eq. m are analogously made invariant using the spurion helds 
E, E and G. 

We remark that the spurion helds in eq. (uni) transform in the same way as the ones for 
the valence-valence four-fermion operators; they diher only in t^r constant hnal values: 
a^Ps ® Pv is replaced by a^Py ® Py in fhe valence-valence case 


D. Chiral Lagrangian 

Assuming that the symmetry in eq. is spontaneously broken down to its vector part, 
the particle spectrum contains light pseudoscalar bosons. These bosons are described by the 

® Spurion fields with the same transformation properties appear also in weak matrix element studies. See 
Ref. 1^ and references therein. 
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field 


E = exp(2!<»//), 


(13) 


which is an element of U{4:N ^ + Ny\Ny). <h is a matrix that collects the psendoscalar helds 
in the usnal way For example, for three sea qnark flavors (m, d and s) and two valence 
flavors (x, y for the valence qnarks and x, y for the valence ghosts) we arrange the fields as 
follows: 


$ = 


/ 

u 

7r~^ 

K+ 

Qux 

Quy 



\ 


7l~ 

D 


Qdx 

Qdy 





K- 

XO 


Qsx 

Qsy 





QL 

-1— 

QL 

X 

P+ 

rL 




Qiy 

Qdy 

Qly 

p- 

Y 

Rly 







Rxx 

Rxy 

X 

P+ 


V 




Ryx 

Ryy 

p- 

Y 

) 


(14) 


Here P~^, X and Y are the xy, xx, and yy valence bonnd states, respectively; P"*", X and Y 
are the analogons combinations of valence ghost qnarks. Rxx is the (fermionic) bonnd state 
XX, and similarly for Rxy, Ryx, and Ryy. Likewise, Qpv represents the mixed bonnd state 
Ph, where P is a sea qnark, P G {u,d,s}, and n is a valence qnark, v G {x,y}7 Qpv is a 
4x1 matrix in taste; while = Q^f is a 1 x 4 matrix. The sea-qnark bonnd state fields 
are U, vr’*', etc. These are 4x4 matrices when we take into acconnt the taste degree of 
freedom. We write 

16 




(16) 


b=l 


(and similarly for 7r~^, ,...) where 

F — {^5; iii.il.. iiilil] 


(16) 


are the sixteen taste matrices in the form of Enclidean gamma matrices (P denotes the 4x4 

Q 

, we inclnde a factor of 2 in eq. (Ild|l and divide Tf, by 2 in 
eq. m in order to keep a consistent normalization of all helds in eq. dn. 


^ We have not bothered to name the mixed bound states of sea quarks with ghost valence quarks in eq. da 
because such states will not enter into the calculations below. 
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Under chiral symmetry transformations in eq. (jHI), the held S transforms as 




( 17 ) 


where L e SU{4:Nf + Nv\Nv)l and R e SU{4:Nf + Nv\Nv)r. 

The chiral Lagrangian is expanded in powers of niq and a^, where rriq stands generically 
for either a sea or a valence quark mass. We adopt a power counting that assumes that the 
size of the chiral symmetry breaking due to the quark masses and the discretization effects 
are of comparable size, he., 


^QCD 


where Rqcd denotes the typical QCD scale, of the order of 300 MeV. A different power 
counting is necessary if one of the two parameters mq/K qcd and is much larger 

than the other one. However, the approximate equalitv in eq. ca is realized in current 
lattice simulations using improved staggered fermions |l 0j |. 

Assuming eq. m, the leading order chiral Lagrangian contains the terms of 0{p‘^, niq, a?) 
and is of the form 


+ MSt) + + a^V. 

o 4 0 


(19) 


Here (...) denotes a supertrace in flavor space and the parameters / and B are undetermined 
low-energy constants.® For our concrete example of three sea quark flavors and two valence 
flavors the dia gon al mass matrix M is given by M =dia.g{mu^j, rrid^j, mg^i, m^, niy, rrix, my). 

As in Ref. |28[ 22, |^, for convenience we leave explicit the mg term, which is allowed 
because of the anomaly, and we take m^ —>■ cx) at the end. Note that 


{<!>) = 2Ui + 2Di + 2Si + X+ Y-X-Y , 


( 20 ) 


where Ui is the taste singlet component of U (eq. (fTH|) L and similarly for Di and Sr 

The potential V in the leading order Lagrangian comprises all terms proportional to oS. 
For our mixed theory it can be written as a sum of three terms. 


V = Us+U'o+Uy. 


( 21 ) 


We adopt a normalization that corresponds to a tree-level pion decay constant / w 131 MeV. 
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The first two terms are just the known taste breaking potentials for the sea sector 271 129||: 


-Us = + h.c. 


( 22 ) 


+ C. 


Y,[{UPs^i^uPs^) + h.c.] 




and 


U's = [(e>sS)(e> 5 S) + h.c.] [(i 5 P 5 S)(e 5 .PsS) + h.c.] (23) 

U V 

+ ^(lPs^}(lPs^^} + C'sA^ ^(UPs^}(LPs^^}- 


Here we introduced the short hand notation ^5 = (^ 5 ,taste ® Iflavor) © ly, etc., for the multi¬ 
flavor generalizations of the taste matrices Tb. For instance, in our concrete example we 
have .^5 = diag(.^ 5 ,.^ 5 ,(^ 5 ,1,1,1,1). The coefficients Ci are low-energy constants. 

Note that all terms in the two potentials Us and U^ only involve the helds in the upper 
left sea-sea block of E. This is easily seen by hrst noting that the matrices Tf, commute 
with the projector P 5 . Consequently, the structure {TbPsETbPsE"^), for instance, can also 
be written as {TbPsEPsTbPsE^Ps). This is not surprising. The taste matrices Tb played the 


role of spurion helds in the derivation of the potential Us +Uq in Ref. 3,|39|. As explained 


in section mn the spurion helds need to be sandwiched by the projector P 5 in the mixed 


theory, he., Tb —>■ PsTbPs- In fact, that was the way we o 
without repeating the details of the spurion analysis in Ref. 


stained the potential Us + U'^ 


m 


39|. 


The remaining potential Uy comprises all terms that stem from the operators in the 
Symanzik ehective action that involve valence helds. Some of these arise from the mixed 4- 
fermion operators in eq. ®. The corresponding chiral operators are constructed by forming 
invariants involving one of the spurion helds in eq. m together with arbitrarily many E 
and Sh Insertions of the mass matrix and derivatives are excluded since they lead to terms 
that are necessarily of higher order in the chiral expansion, at least 0{p^a^^mqO?). 

We can only form non-trivial invariants with the helds P and G, which collapse to the 
same term once the spurion helds are set to their hnal value: 


(PiSP2St)^a2(r3Sr3St), 
(G1SIG2S) ^ (rsSVsS) . 


(24) 
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Here we used P 5 = ^(J + T 3 ) and Py = \{,I — T 3 ), with = diag(/ 5 , —ly) and dropped 
an irrelevant factor of 1/4. When we write (J ± r^) for Pi 2 and the helds S and St 
cancel whenever they sandwich the identity matrix. The only non-trivial operator is the one 
involving two matrices. 

The last terms in Uy stem from the four-fermion operators involving only valence helds.® 
As we remarked at the end of the last section, the corresponding spurion helds transform 
exactly as the ones for the mixed four-fermion operators. The only diherence is the hnal 
constant value; a?Ps ® Py is replaced by a^Py ® Py. This change only involves a sign hip 
in the hrst projection operator, and therefore leads to the same term (rsErsEl) for the 
chiral Lagrangian.^® 

We conclude that the leading order chiral Lagrangian for the mixed action theory with 
staggered sea and GW valence quarks contains only one more operator compared to the 
chiral Lagrangian of SXPT: 

Wv = -GMix(r3Sr3St). (25) 

The potential V for the mixed action theory involves nine unknown low-energy constants 
compared to eight in SXPT. We want to emphasize that the eight constants in Us + Wg 
are the same constants as those in SXPT. Some combinations of them have already been 
determined by hts to staggered lattice data [k| . 

The potential Us +Wg breaks the S'P(4) taste symmetry in the sea sector but not entirely 
- an accidental S'0(4) subgroup remains unbroken [^. The part Uy, on the other hand, 
preserves the full SU (4) taste symmetry. This is expected because the four-fermion operators 
in eq. which give rise to Uy, are trivial in taste space. Interaction vertices involving 
pseudoscalars with one or more valence-quark constituents stem from Uy only. The SU (4) 
taste symmetry implies that correlation functions that include such external pseudoscalars 
respect the SU^A) taste symmetry in one-loop NLO diagrams. Analytic, taste symmetry 


® The allowed valence bilinears either have the same symmetries as the lowest order terms, in which case 
they merely give 0{a^) corrections to lowest order parameters, or they have different sym metries (violate 
0(4) rotation invariance) and contribute only at higher order in the chiral expansion |l8l| . 

This also explains that this term is present independently of a change of basis in the Symanzik effective 
action. Using eq. one could replace either the mixed or the pure valence four-fermion operators at the 
expense of (4>r'I')^, which is invariant under transformations in Gpqqcd- Nevertheless, one can remove 
only one type of four-fermion operators, the other type still gives rise to the term (rsETaSl) in the 
chiral Lagrangian. 
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violating contributions do appear at NLO: a^niq)-, while non-analytic symmetry 

violating contributions start beyond NLO. 


III. PSEUDOSCALAR MASSES AND DECAY CONSTANTS 


In this section we compute the one-loop expressions for the (valence-valence) pseudoscalar 
masses and decay constants. For concreteness we now restrict ourselves to the case with three 
sea quark flavors and two valence quark flavors. This is the most relevant case phenomeno¬ 
logically. Further, the resulting one-loop expressions can be readily used in the analysis 
of unquenched conhgurations generated with the Asqtad action by the MILC collaboration 
(see Refs. j^[l^ and references therein). 

In the following we also adjust for the so-called “fourth root trick”, which is commonly 
employed in staggered simulations in order to reduce the taste degree of freedom from four 
to one. In the context of the chiral effective theory this adjustment requires proper insertions 
of factors of 1/4 in the sea quark loop contributions in our expressions [2^, depending on 
the quark flow that corresponds to the meson loop diagram in the chiral effective theory. 

The fourth root trick raises legitimate locality questions and its validity is controversial. 


Recently various s 


ods 


udies have addressed this issue using either numerical or analytical meth- 
In the following we assume that the “fourth root trick” 


3,Q,Q,y,Q,Q,0,Q 


can be given a he 
described in Ref. 


d theoretically sound underpinning, so that we can follow the procedure 

m. 


A. Leading-order masses and propagators 

At tree level, the new operator of the mixed theory, Uy, contributes only to the masses of 
valence-sea mesons, represented in eq. (HI by Qfv (-^ ^ n G {x,y,x,y}). Sea-sea 

and valence-valence mesons get no such contributions because a block-diagonal $ commutes 
with Ts, and SSl = I. Similarly, the potentials in the sea sector. Us and eqs. ()22II23|1 . 
give no contribution to the valence-valence mesons. Thus the valence-valence mesons obey 
the continuum-like mass relations exemplihed by 

rup = B{mx + my) . ( 26 ) 
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Of course such relations follow immediately from the exact chiral symmetry (for massless 
quarks) in the valence sector. 

In the sea-sea sector, lAs and Wg contribute, and the tree-level results are identical to 
those in Ref. For a meson of taste b made up of sea quarks F and F' [F ^ F'), we 
have 


m 


= Bimp + m'p) + a^A(^b), 


(27) 


with 


A(e5) = Ap = 0, (28) 

A('Cm 5) = Ayi = — {Cl F 3O3 -1-04-1- SOg) , 

A{i,,) = At = ^ (2O3 + 2O4 + dOg) , 

A(^^) = Ay = — (Oi -|- O3 -|- 3O4 -|- 30 g) , 

A(6) = A, = ^(403 + 404 ). 

As mentioned above, Uy contributes to the masses of the valence-sea mesons. For the 
Fx meson, with held Qfxi the mass is 

= B{mF + m„) -h o^Amix , (29) 

where 

Am,. ^ . (30) 

The violation of the Goldstone theorem in eq. (El clearly arises because there is no lattice 
axial symmetry that mixes valence and sea quarks. 

In a simulation using staggered sea quark conhgurations, Amk could be directly deter¬ 
mined from the propagator of a mixed meson with one GW valence quark and one staggered 
valence quark.Such a direct determination of Amix would be useful because it would 
reduce the number of free parameters in chiral-log hts. For example, Amix enters (through 
into the NLO expression for the decay constant of a meson with two GW valence 
quarks (see section mini. 


At least within the context of staggered chiral perturbation theory, taking the fourth root of the staggered 
sea-quark determinant will not change the standard equivalence between masses of particles on internal 
and external lines. The fourth-root procedure only changes the weighting, not the mass, of an internal 
meson made from one valence and one sea quark. 
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Since we have no a priori reason to expect a particular sign for Cjuix (or equivalently 
Amix), eq. (j2ni) shows there is a possibility of a lattice “Aoki phase” if Amix < 0. This 
would be similar to the type of lattice phases for staggered quarks discussed in Refs, 

1^ . The direct measurement of Amix discussed in the previous paragraph would determine 
if this possibility is realized in practice. 

The flavor-charged (non-diagonal) helds in eq. ()14|1 have only connected propagators in 
the quark-flow sense, Fig. ^ while flavor neutral (diagonal, e.g., U or X) mesons also have 
disconnected contributions. Fig. |21 The only complication in the connected case is getting 
the sign of the propagator right for mesons with one or more ghost valence quarks. Using 
[AB]{p‘^) to denote the Euclidean space propagator of held A and B with momentum p, 
examples of connected propagators are 

1 


[P+F-] (p 2 ) = 
(p^) = 


[XX] 

XX 


QuixQ 




'^y ^xy 


(/) = 
(/) = 

conn 

(p^) = 

(/) = 


9 I 9 ’ 

p^ + rrip 
Sab 

p2 _j_ ’ 

1 

9 i ^ ’ 

P^ + TTlx 
-1 

9 i ^ ’ 

p^ + m% 


(31) 


t 

UjX 


s., 




pi 




(p^) = 


P + ^xy 

Here a,b = {1,2 ,..., 16} are meson taste indices as in eq. (ITH|l : while i,j = (1, 2, 3,4} are 
quark taste indices. 


< 


> 


FIG. 1: Quark flow diagram for a connected meson propagator. 


Disconnected propagators for havor neutral mesons can be generated by the anomaly (tUq) 
term in eq. (HI, which gives a “hairpin” interaction (two-meson vertex with disconnected 


quark how). Summing such graphs as in Refs. 




gives, for the disconnected X-Y 
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+ < , * > 

+ < ^ " - > 

+ * * • 


FIG. 2: Quark flow diagrams for the disconnected meson propagator. 


propagator 

[^>^]disc ip") 


_ jP^ + + ^Di){P^ + ^%) _ 

3 {p^ + m\){p^ + mY){p^ + 


(32) 


where, for concreteness, we have again assumed three sea-quark flavors. The vr®, pj and t]j 
are the mass eigenstates in the flavor-neutral, taste-singlet channel, found by diagonalizing 
the mass matrix including the rUg term. The subscript “disc” is included for clarity, but 
of course the XY propagator has no connected contribution. In the sea-sea sector, only 
the taste-singlet, flavor neutral mesons feel the anomaly hairpin vertex. In addition, sea-sea 
flavor neutral mesons of vector or axial taste get hairpin contributions from j29| . Since 
the corresponding disconnected propagators do not enter into the quantities calculated in 
sections nm] and nun we do not write them explicitly here. 

It is convenient to take ttIq — > cx) at this point to decouple the r]'j. In the case of most 
current interest, with three staggered flavors and the fourth root of the determinant taken to 
eliminate the extraneous taste degree of freedom, for large mg Q. Equation 

then becomes 


[^>^]disc ip") = -F 




3 + mY){jp + rn?+ m. 


(33) 


7 ) 1 / 


Other disconnected valence-valence propagators are found from eq. o by substitution: Let 
Y —7 X for X ^ Y for If = rrid, as in the MILC simulations P. [ic|. 
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then 



( 34 ) 


If we take the fourth root but keep the number Nj of flavors arbitrary, then 
A^/tUq/S, and the more general version of eq. ()dd|l is 



(35) 


where L runs over diagonal flavor neutral states {U, D, S, ...), and L' run over the neutral 
mass eigenstates (tto, t], ...), excluding the rj'. 

B. NLO valence-valence mass 

We are interested in computing the one-loop correction to the meson made from valence 
quark x and y, i.e., a P~^. The self energy comes from the meson graphs in Fig. El We 
thus need the four-meson vertices generated by the kinetic energy, by the mass term, and 
by Uv, with at least one P~^ and one P~ held. 



(a) 


(b) 


FIG. 3: Meson graphs for the P'^ self energy. Graph (a) has a connected internal propagator 
(Fig. n in the quark flow picture); while graph (b) has a disconnected internal propagator (Fig. Ej)- 

Expanding the kinetic energy, and including the minus sign for a vertex, gives the follow¬ 
ing terms with derivatives acting on both and P~: 



(36) 


F 
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The fields in this expression are defined in eq. m- The summation index F runs over 
sea quarks, typically u, d, and s. There is also an implied sum over the taste index in the 
product Q^Q. As usual, terms with a derivative on or P~, but not both, are not needed 
since the corresponding diagrams in Fig. |21will vanish by symmetric momentum integration. 
The kinetic energy terms with no derivatives on the P^, P~ fields are 


p^(2) _ 1 p+ p- 

^KE - T72 ^ ^ 




\d,X)‘‘ + (d^yy + + d,Rl d,Rf, + d^Rl, d,R,, (37) 

df^Ryy + ^ (SpOk + SpOfi k'SFj) - 'id^x d^Y 


Similarly, the mass term and the “mixed potential” Uy give, respectively. 




B 

3f 


p+p- 


{rrix + my)P^ P + (Sma, + my)X‘^ + {mx + 3my)Y‘^ 


(38) 


T(3?7ra; T TTly'jP^^Pxx T 2 {jTlx p A ^xy^^y^ A ip^x P ^'^y')^yy^y' 

P ^ |(2ma; PniyP mp) Q^^Qfx + P 2my P mp) Q^yQpy^ 


p2{mx P my)XY 


and 

Vmix = P*P- {Qf.Qf, + Q'fvQf^) . ( 39 ) 

where we have used eq. (PU|). 

By including the ghost quark contributions we have guaranteed that quark loop terms 
from the valence quarks will be canceled automatically, thereby accomplishing the partial 

Q nching. However, we still must understand the meson diagrams at the quark-flow level 
in order to adjust for the effects of the fourth root procedure on the staggered quarks. 
If we assume three flavors of sea quarks for definiteness, the fields P^ and P~ correspond 
to ^45 and ^ 54 , respectively ^cf. eq. (|14|) j. The vertices, Eqs. (j3bj) , (j37j) , (j3Sj) and (Ei, then 
come generically from two types of terms in the supertrace of four d) fields: terms where 
P^ and P~ are next to each other i Y,i (or ~ X]. <F 54 $ 45 <I) 5 i$i 5 ), and terms 

where and P~ are separated: ~ $ 45 $ 55 $ 54 $ 44 . These can be represented by the vertices 
(a) and (b) in Fig. Irrespectively. The index i in Fig.|31(a) should be summed over all valence, 
ghost valence and sea quarks. The last contributions (proportional to XY) in Eqs. 

(E7|> and dSHD correspond to vertex Fig. Hb); while all other terms come from the sum over 
i in Fig. EKa). 
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Note that the graphs in Fig. H] represent quark flow only; the numerical value of each 
graph depends on the term in the Lagrangian that generates it, and may also depend on 
the free index i. In particular the mixed potential Uy generates only vertex Fig. EJ^a) terms, 
and only gives non-zero coefficients of terms where i is a sea quark flavor. This follows from 
the fact that in eq. (ESI) is proportional to the identity in the pure valence sector, and 
therefore lAy reduces to a field-independent constant in that sector. 



X 


y 



(a) 


(b) 


FIG. 4: Quark flow graphs corresponding to the four meson vertices in the the P'^ self energy, 
Eqs. (EH), EH: (EH and EH- The horizontal x, y lines produce the external fields. Graph 
(a) represents terms where the P~^ and P~ are next to each other in the supertrace; an almost 
identical graph with x <-> y is not shown. The free index i represents any quark type, but the 
numerical coefficient of the graph may depend on i. Graph (b) represents terms where P~^ and P~ 
are not next to each other in the supertrace. 


When two meson lines at the vertices are contracted, as in Fig. El the quark flow diagrams 
in Fig.Elresult. The connected contraction of Fig.EI^a) gives Fig.l^a), which clearly involves 
an internal sea quark loop. This means that terms from Fig. I^a) where i is a valence or 
ghost valence quark must cancel in the connected contractions. This arises algebraically 
from Eqs. EH, EH and using eq. EH): Contractions of B)R terms cancel connected 
contractions of XX and YY and P^P~ terms. This leaves just connected Q^Q contrac¬ 
tions (Fig. Efa)), disconnected XX and YY contractions (Fig. Efb)), and disconnected XY 
contractions (Fig. Efc)). 

At this point it is easy to make the “by-hand” adjustment necessary to correspond with 
the fourth-root procedure. The only explicit sea quark loop in Fig. Elis in diagram (a), so we 
just divide those terms by 4. We already know such terms come only from contractions; 
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(a) (b) 



(c) 


FIG. 5: Quark flow graphs corresponding to the P~^ self energy, Fig. |3J Graphs (a) and (b) come 
from connected and disconnected contractions, respectively, of the internal meson lines in vertex 
Fig. IHa) ; graph (c), from the disconnected contraction in vertex Fig. Hb). Iterations of sea quark 
loops in the disconnected propagators, as in Fig. [2 is implied in graphs (b) and (c). 


dividing them by 4 is equivalent to ignoring the implicit sum over the 4 tastes in the Q^Q 
terms. There are also sea quark loop insertions in the disconnected meson propagators in 
Fig.El^b) and (c). Such insertions have already been corrected for the fourth root in eq. ( 1221 ) 
or eq. 

Let Ep+ be the self energy, dehned to be the negative of the sum of self energy 
diagrams. At NLO, we have 


Sp+(/) 


1 


((To rup + (Ti + analytic terms , 


(40) 


where (Tq and ai come from the one-loop diagrams and are independent of momentum p. 
Since we have not determined the NLO chiral Lagrangian, we cannot express the analytic 
contributions in terms of chiral parameters. For the quantities of interest, however, it will 
be sufficient for our purposes to write down the most general contributions consistent with 
the symmetries. Putting together the one-loop contributions from the vertices Eqs. dSHD, 


21 













and (jnni), we get 


o-Q = - 


E 


1 1 
+ 


+ rrip^ + m 


( 41 ) 


Fy , 


+ [XXl^^iq^) + [YYl,Jq^) + 4 


CTi = - 


E 


1 1 
+ 


^ + q^ + mj,y 

+ [^■>fl<ll.c(9") + - 2 |2fyld,.c(9") 

where, as usual, F runs over sea quark flavor {e.g., u,d,s), and 

f = IGtt^ 


(42) 


(43) 


( 2 ^)* ’ 

suitably regulated. The disconnected propagator is given in eq. dMl) or eq. (jsni)- In eq. (01, 
we have dropped (constant) quartic divergences and have used Eqs. (^1 and P7jl to 

replace factors of quark masses with the corresponding meson masses. The simple identity 


+ ^ Jdiscl^ ) ^2 ^^2 




(44) 


q^ + 

has allowed us to remove all explicit factors of q^ in the integrands involving disconnected 
propagators. 

We now focus explicitly on computing the mass. The chiral symmetry in the valence 
sector implies that the analytic contributions to (mp^°)^ must be proportional to the tree- 
level rup oc {rrix -|- Fiy). In the continuum limit, these contributions just go over to the 
standard form in terms of the Gasser-Leutwyler parameters Li Q. At lattice™ 
the only possible new NLO analytic term is Ca^m?p, where C is an unknown constant that 
depends on the details of the lattice action in both the valence and the sea sectors. Indeed 
it is easy to write down terms that will appear in the NLO Lagrangian and will contribute 
to C, for example. 


a-^ {PvdyY.PydyL^) , ( 45 ) 

a^Us{dyY.dyL^) , 
a^Us {T.M^ + MT)) , 
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with lAs given in eq. (j22I)- The one-loop corrections to the mass squared are then found 
by evaluating Sp+(p^) at p^ = —m?p, giving 


[m 


NLO\2 


B{m^ + rriy) 


= 1 + 


(c^o - cTi) + ^ (2L8 - Ls) (m^ -1- niy) 

^ 32B 


(46) 


{ 2 Lq — L 4 ) ^ ^ 771 j- P a^C, 


where, from Eqs. (14211 and 


ao - = -2 / [XY],^^{q^) . 


(47) 


The integral in eq. (14711 c an be evaluated in terms of the chiral logarithm and residue 
functions dehned in Refs. (291 130 1. For completeness we include the dehnitions here. The chi¬ 
ral logarithm functions, coming from integration over a single and double pole, respectively. 


are: 


i{rn^) = rn^ ( In ^ + Si{mL) ) , 
e{m^) = - ( In ^ + 1 I + ( 5 s(mL) , 


(48) 


where is the chiral scale, L is the spatial dimension, and the hnite volume correction 


terms are 


5i{rnL) = 4^ 


r^O 


i7i(|r|mL) 
7 nL\f\ 


(49) 


63 (mL) = 2 '^ Ko(\f\mL) . 

r^O 

Kq and Ki are Bessel functions of imaginary argument, and r, which labels the various 
periodic images, is a three-dimensional vector with integer components. 

The residue functions allow one to write integrals over ratios of products of {q‘^+m‘^) 
terms, such as [-AF]jjg^(g^), as integrals over single poles. They are dehned by 






(50) 


The residues are a function of two sets of masses, the “denominator” set {Ad} = 
(mi, m 2 ,..., mn} and the “numerator” set (p) = (pi,/i 2 ,...,/ifc}- The indices j and z. 
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1 < i, * ^ '^5 refer to particular denominator masses; the prime on the product in the de¬ 
nominator of eq. CT means that i = j is omitted. In cases of degeneracy, we also need the 
double-pole residue functions, 


^ ■ ( 51 ) 

We can now write out eq. (m explicitly in various useful cases. In the Nf = 3 partially 
quenched case, with no mass degeneracies, we have 


d 


,[n,k] 


(m 


NLO'i2 
P+ ) 


B {rrix + my) 


= 1 + 






(52) 


i=i 


165 , , , , 325 , , , ^ 2^ 

H—^ (2L8 - Ls) [mx + my) + —p^ [2 Lq - L 4 ) {mu + md + ms) + a C. 


The index j is summed over the denominator masses, as it will also be in subsequent cases; 
the mass-set arguments are 

= {mx,mY,m^o,my^} , 

= {muj,mDi,msj} . 


(53) 


Equation (1221), is identical to the corresponding continuum partially quenched result with 
no degeneracies [^, except for the explicit discretization term a^C and the fact that the 
neutral sea-sea mesons are specihed to be taste singlets. For us, “no degeneracies” means 
that none of the seven meson masses in eq. (El are equal. We have chosen the normalization 
of quark masses in the chiral Lagrangian, eq. dH, so that the same constant 5 appears 
in the relation between quark and meson masses in both sea and valences cases. However, 
degeneracies between sea-sea and valence-valence mesons in eq. (El would not imply equal 
valence and sea quark masses, because of the splitting of the taste-singlet sea-sea mesons, 
eq. (j21. 

Taking = m ^ muj = moj = m^o, but with no other degeneracies, then gives 


[m 


NLO'i2 
p+ ) 


5 {mx + my) 


= 1 + 






(54) 


i=i 


165 325 N / \ 

{ 2 Ls - L5) {mx + my) ^ {2Lq - L4) (2m -f m^) a C , 




with 


{M 




xy,/} = {'mx,mY,my^} , 
{nf} = {muj,ms,} . 


(55) 
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When niu = mj, = m and tiIx = my mx = my, bnt no degeneracies between sea-sea 
and valence-valence mesons, we have a “partially quenched pion” with 


(m 


NLO'i2 
p+ ) 


2 Bm^ 


= 1 




[ 2 ] 


,| 21 - 


(56) 


+ DM({>fg,};{db) «(>"?) + 

J =1 / 

(2-^8 ~ -bs) H— -p^ {‘^Lq — L 4 ) ( 2 m -\- m^) + a^C , 


/ 


P 


where 


= {mx,myj} , (57) 

[91 BBBV 

and {/i/ } is given by eq. 

At hnite lattice spacing, the cases that most resemble the full (unquenched) theory are 
ones with degeneracies among the valence-valence and sea-sea mesons in eq. diSI)- For 
current purposes, we might call a “full pion” one with mx = my = mjj^ = my,^ = m^o, 
which requires = my and mu = m^ = rh, but mx > m, since the taste splitting A^o in 
eq. (| 2 ID is positive.In this case, we have 


(^NLO)2 

2 Bm^ 


= 1 + 




( 68 ) 


(SLg — L,) H—p,- {2Lq — Li) (2m + mi + a^C , 


I J .2 ^ 0 / v“"-x/ I J .2 

where we have used eq. dSl to simplify the residues. Similarly, calling a “full kaon” a meson 
with mx = mjjj = moj and my = msj, results in 


(^NLO)2 


B{mx + my 


(2Ls - L 5 ) (mx + my) -F 


(59) 


32B 


I J .2 v“^o ^0/ \"-x 1 '—y/ 1 J .2 

Equations dSHl) and (jsni) clearly approach the standard results 


{2Lq — L 4 ) (2m -|- mg) + PC . 
as a —>• 0. 


C. NLO valence-valence pseudoscalar decay constant 

The decay constant fp of the meson is defined by the matrix element of the corre¬ 
sponding axial current, 


(0 \jy5\Pip)) = -ifpPi ■ 

This is at least true for the simulations carried out by the MILC collaboration. 


( 60 ) 
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In the LO chiral theory of eq. (P|) . the axial current is given 


3 lib 






( 61 ) 


Here A projects out the appropriate flavors: With three sea-quark flavors as in eq. (IT^ . the 
valence quarks x and y correspond to indices 4 and 5 of d*, and then A^ = 

At NLO, the decay constant has the form 


fp 


NLO 


= 1 + 


;6fp + analytic terms . 


( 62 ) 


/ ■ ' leirV" 

The term 6fp comes from the one-loop diagrams, and the analytic contributions are gener¬ 
ated at tree-level by NLO terms in the chiral Lagrangian and corresponding corrections to 
the current in eq. (jnn). There are two contributions to 6fp, 


Sfp = ^ ^ 


( 63 ) 


where comes from diagrams generated directly by expanding the current in eq. p[|) 

to cubic order in 4), and 6Zp is the one-loop wave function renormalization. From eq. iOl) 
we have 

SZp ^ -lePf . ( 64 ) 


As in Ref. 
we then have 


Q, 


dp^ 


a straightforward calculation shows that 5/' 


current 


= —2 6Zp. From eq. la, 


« 3 1 

= r' = “2 


E 


+ 


+ rrip^ rripy 

+ I^A']„„(g^) + - 2 |A'y]„.„(5^) 


(66) 


with F summed over sea quark flavors. 

The analytic terms in eq. (El come only from NLO terms in the chiral Lagrangian 
with derivatives, which affect the decay constant either directly, through wave function 


We use the Noether current corresponding to axial vector rotations as our partially conserved axial vector 
current. This is justified if the analogous current is used in the underlying lattice theory. A convenient 
method to construct this current in numerical simulations is described in Ref. . Alternatively, one may 
employ the corresponding covariant pseudoscalar density to define the decay constant On the other 
hand, using a conserved but not covariant axial vector current would result in extra terms proportional 
to arrix and aruy, which are not captured in the following results. 
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renormalization, or indirectly, by leading to higher corrections to the axial current. Thus, 
O(a^) corrections to the chiral Lagrangian have no effect on eq. (Ib2j) . There will however be 
analytic terms from 0{p‘^), 0{mp‘^), and 0{a^p^) Lagrangian corrections. The former two 
are identical to those in the continuum, and produce terms proportional to quark masses. 
The effects of the latter corrections on fp can be absorbed into a single term proportional 
to a^. We thus have. 


fNLO 

Ip 

f 


= 1 + 


ojp + (m„ + my) + } ^mp + a T , 

J F 




P 


( 66 ) 


where L 4 and L 5 are standard j 2 ^, JF is a new constant, and 6fp is given by eq. (inni)- 
Because Lagrangian terms of O(ma^) do not affect JF, it is easy to see from the discussion 
surrounding eq. (iiD that JF is independent of the corresponding constant C occurring in 
the expression for the meson mass, eq. (j46|l . Like C, T depends on the details of the lattice 
action in both the sea and valence sectors. 

We can now write out the NLO expression for the decay constant in various special cases. 
In the Nf = 3 partially quenched case, with no mass degeneracies, we have 


I 


NLO 

P 


/ 


= 1 + 


167r2/^ 


1 




(67) 


+- I {Ppmmj,) + {Ppmmi.) 


,[ 3 ] 


,(3) 


[ 3 . 3 ], 


,[ 3 ] 


,(3)l 




4 


i=i 


-2 


8B r . , 165 ^ ^ 2 ^ 

L 5 {m^ + my) + } ^mp + a R , 


where F runs over u, d, and s. The mass sets {■M'^xyi) given by eq. 

and 


,(3)’ 


= {mx,m^o,myP , 


( 68 ) 
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With rUu = md = m ^ mu^ = moj = m^o, but no other degeneracies, the result is 


fp 


NLO 


/ 


= 1 + 




[2^ {Kx) + ^ {^Ix) + 2^ {mly) + £ {m%)] 


( 69 ) 




,[ 2 ] 


,( 2 ) 


[ 2 , 2 ], 


,[ 2 ] 


,(2)l 


6 


i=i i=i 

3 

[3-2]aA^[3] 1. r„(2)lA^^_2 


m 




j=i 

85, , , 165, , 2^ 

-—Lg (m^; + my) + -^-^4 (2m + mj + a 5 . 


Here, {■M'^xyi} are given in eq. (jKK|l : while {M.^xi) dehned in eq. (jFfjl (h 

{M^y]i} take X^Y). 

For a “partially quenched pion” with m^, = m^, there is considerable simplihca- 
tion because the disconnected contributions in eq. (jnsi) will cancel. Taking in addition 
rriu = TYid = for simplicity, we have 

^NLO ^ 1 


,( 2 ) 


[ 2 ] 


or 


/ 


= 1 + 


2 ^ (rnL) - ^ ("^y 


167r2/2 

85, , , 165, , ^ , 2^ 

+—Lg (2ma.) + (2m + mj + . 


( 70 ) 


There is no obviously preferred way here to dehne a “full pion” to make the NLO corrections 
take on a continuum-like form. In the a —> 0 limit, the splitting a^Aj^in eq. will vanish, 

and the logarithm terms will clearly approach the standard form 2^: —2^ (m^) — ^ (m^). 
At hnite lattice spacing, we can choose rux so that m^^, and m^^, have the masses of the 
sea-sea pion and sea-sea kaon of any one taste, but there seems to be no advantage in doing 
that. In particular, the value of m^. so chosen will not be the same in general as the value 
needed to give the logarithms in the meson mass their continuum-like form, eq. 


For the kaon, it makes some sense to dehne a “full kaon” as we did in Sec. unBi mx = 
muj = moj = m^o and my = msj. This at least gives the disconnected contributions the 















form they would have in the continuum. We then have 



In the continuum limit, niux = '^^0 = rUj, msx = 
thereby reproducing the known result 22 j |. 


rriK, TTtrjj = and rUsy = msj, 


IV. DISCUSSION 


Our results for 2+1 sea quark flavors are currently the most relevant ones, since they 
can be applied to simulations using the existing conhgurations generated by the MILC 
collaboration. Eqs. (Hi and (Hi describe the quark mass and lattice spacing a dependence 
of the pion masses and decay constants, and these expressions can be directly htted to lattice 
data obtained with Ginsparg-Wilson valence fermions. 

The number of unknown £t parameters in these expressions is fairly small. For instance, 
the pion mass depends on the usual low-energy constants of continuum XPT (/, B and two 
familiar combinations of Gasser-Leutwyler coefEcients), the sea-sea meson masses 
the sea quark masses rh, rUg and the constant C. At one loop order we can express the sea 
quark mass combination 2m through leading order sea-sea meson masses (cf. eq. (I27|) h 

B{2m + ms) = ^m^++ m^+. (72) 

These masses as well as m^o and m^^ have already been measured by the MILG collaboration 
and are therefore not unknown parameters. The only true unknown parameter in addition 
to the ones from continuum XPT is thus the constant C. 

The measurement of the singlet meson masses is difficult because disconnected diagrams contribute to 
the correlator. For = rrid, however, there are no disconnected contributions to the ttJ propagator, and 
its mass is degenerate with the tt/ mass. The rji mass is also not affected by disconnected diagrams if 
m„ = my + rus and the limit mo —> oo is taken. In that case it is consistent to employ eq. m where 
the ms, and mu, masses are from connected diagrams only. In reality m§ is not infinity, and there can 
he rj — r]' mixing, which would be proportional to (m^ — mu)lmQ. Such corrections are not taken into 
account in the MILC determination of m^,. 
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Similar statements apply to the pion decay constant in eq. dZOD. Even though the masses 
of the valence-sea mesons have not been measured yet, they can be straightfor¬ 
wardly determined from the propagator of the mixed meson and a linear fit to the leading 
order mass formula in eq. (El- Using this information leaves one additional parameter, the 
constant JF, besides the familiar parameters of continuum XPT. 

Let us compare our results for Ginsparg-Wilson valence quarks with the corresponding 
expressions for staggered valence quarks. The one l oop expression for the Goldstone pion 
TT^ in the 2-1-1 flavor case reads (see eq. (75) in Ref. j29|) 

,NLO\2 
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The first two rows of this expression give the corresponding result in eq. dSHD for Ginsparg- 
Wilson valence quarks. The remaining contributions involve many more sea-sea meson 
masses as well as the two “hairpin” parameters 6 y and These parameters cannot be 
expressed in terms of leading order charged meson masses and are therefore true uncon¬ 
strained fit parameters. 

The one-loop result for the pion decay constant - given in eq. (27) in Ref. Q - also has 


These two parameters are combinations of the low-energy constants in the potential W: 5y = C^v — Bsv, 
and analogously for [ 2 ^ . 
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contributions proportional to Sy and 
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The factor of 1/16 in the hrst line is canceled in the continuum limit by the sum over b, 
which runs over all sixteen different meson tastes. As was the case for the pseudoscalar 
masses, the corresponding expression for the mixed theory, eq. (iruj). is much simpler and 
does not involve the contributions proportional to the hairpin parameters 6 y and 5 ( 4 . Note 
that the constants C, T in eqs. dzg) and dZlD are different from C and T in eqs. (P|) and 

(|7n|). 

Obviously the functional dependence in the expressions for staggered valence quarks is 
much more complicated and involves more undetermined parameters than in the corre¬ 
sponding results for the mixed theory. However, the physical low-energy constants, the 
Gasser-Leutwyler coefficients, enter the expressions in the same way. Therefore, as already 


pointed out in Ref. ^|, mixed simulations may be used to extract these phenomenologically 
relevant parameters from numerical lattice simulations. 

The NLO formulas computed in Secs. 1111 Ri and lTTl Gl also make concrete a rather obvious 
fact about mixed action theories: At non-zero lattice spacing there is no way to dehne 
equality of valence and sea quark masses in order to have all properties that might be 
desired of a “full” (unquenched) theory. The lattice theory will always have some features 
of partial quenching, and it is only in the continuum limit that the full theory is obtained. 
Depending on the purpose one may wish to choose various dehnitions to match the sea and 
valence quark masses. Since the scalar correlator is very sensitive to the effects of partial 


quenching 


53j |. it has been proposed to choose the valence quark mass such that the scalar 


correlator does not have a negative contribution 


For the results reported in Ref. 


Q 
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the masses were chosen so that the valence-valence pion mass coincides with the Goldstone 
pion mass made of staggered quarks, he., = w? + . The NLO results have suggested yet 

TTg 

another dehnition with m^+ = '^^+- In this case, as we have seen in Sec. mm some chiral 
logarithms resemble their continuum form and one might expect smaller partial quenching 
effects than with other dehnitions.^® From a theoretical point of view all these dehnitions are 
equally good since they guarantee that full unquenched QCD is reached in the continuum 
limit. Practically, they differ with respect to the size of the partial quenching effects at non¬ 
zero lattice spacing, and the quark mass tuning can be rather difficult to achieve, depending 
on, for example, the statistical errors in the observables used for the matching. 

Nevertheless, at least in the context of chiral perturbation theory, there is no fundamental 
difficulty with using mixed action theories to simulate QCD. The effects of hnite lattice 
spacing can be controlled by hrst htting to the chiral forms of the type derived here, and 
then extrapolating to the continuum limit. Furthermore, as in the pure staggered case j^, 
we expect the chiral and continuum limits to commute in the mixed theory for any quantity 
that has a well-dehned chiral limit in the continuum. 
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With the definition in Ref. the scalar correlator becomes negative [^, a clear signal for partial 
quenching. 
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APPENDIX A: CONSTRUCTION OF MIXED FOUR-FERMION OPERATORS 


In this appendix we construct the mixed four-fermion operators in eq. © that enter the 
Symanzik effective theory at We closely follow the procedure and notation in Ref. 271 

where the four-fermion operators made of sea quarks only were constructed. The method 
determines hrst all lattice four-fermion operators without derivatives and mass insertions 
that are singlets under the symmetries of the lattice theory. Taking the continuum limit of 
these terms results in the allowed continuum operators that appear in the Symanzik effective 
action. 

First we convert the staggered helds into hypercube helds, since these helds yield the 
proper continuum helds when the lattice spacing is sent to zero. Following Refs. ^ 
the lattice is divided into hypercubes containing 16 sites whose coordinates are written as^^ 


= 2y^ + T]^. (Al) 

The hypercube vector r] labels the sites within the hypercube and its components are 
either 0 or 1. In terms of the site variables y, x ^-nd the gauge links U the hypercube helds 
are dehned by 

i^s,aa{y) = ^'^'^T^i‘^y,‘^y + v)xs{‘2y + v), (A2) 

V 

^s,aaiy) = ^5^Xs(2|/ + h)f^^(2|/,2r/ + r/)F;"“, 

V 

where U{2y, 2y + ri) denotes a product of link variables along a path going from 2y to 2y + ri, 
and 


r. = • (A3) 

The indices a and a represent the Dirac and taste index, respectively (we suppress the havor 
and the color index). 

Using the hypercube helds we now construct all mixed four-fermion operators 0/^i that are 
allowed by the symmetries and that correspond to dimension 6 operators in the continuum 
limit. Four-fermion operators that contain derivatives and/or quark masses are higher than 

We use lattice units and set a = 1 in this appendix. 
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0 (0^) in the Symanzik action and can be ignored here. Since no quark mass appears the 
operators must be invariant under the full chiral symmetries of the massless lattice theory. 
The construction proceeds in hve steps 0: 


1. Multiply a sea quark bilinear on a hypercube by a valence bilinear at the same lattice 
point and sum over all hypercubes, 

04f(r5,ry) = (A4) 

y 

The flavor symmetries dictate that the sea quark bilinear is an SU{Nf) singlet while 
the valence bilinear is a singlet under SU{Nv\Nv)- T^ represents an arbitrary tensor 
product 7 a ® ® of a gamma matrix acting in Dirac space, a color generator and 

an SU{4:) taste group generator Similarly, Ty denotes an arbitrary combination 
'jB <8) acting on the valence helds. It does not include a slot for a taste matrix since 
the valence helds do not have the taste degree of freedom. 

In eq. flA4jl all the sea quark indices are contracted with T s while all the valence indices 
are contracted with Ty, so the operator truly is a product of two bilinears. One can 
write down other operators that do not have this simple structure. For example, one 
could contract the color indices of 'ijjs and tfjy and the indices of V’v and '0s- Similarly 
one can contract the Dirac indices in such a “twisted” manner. However, all these 
ope.at„. a.e .edaadant Q. MaU.ng use of Fie. .deut.ties one can always -untwisf 
these operators and bring them into the form in eq. ()A4fl . 

On each hypercube there are 16^ possibilities to form a valence held bilinear. The one 
chosen in eq. (IA4fl involves only the helds at the lattice point y where the staggered 
hypercube held lives. Close to the continuum, all other valence bilinears can be written 
as the one in eq. ()A4fl . plus terms involving derivatives, which we can drop. 


2. The sum over all hypercubes in eq. ()A4jl makes the operator 04 f invariant under lattice 
translations by one hypercube, i.e., y ^ y + 1. In order to obtain the part that is 
invariant under single site translations we apply the projection operator 

-p = (A5) 

y 

with Tg (g) Ty being the translation operator in the /i-direction. The translation oper¬ 
ator acts diherently in the sea and valence sector. In the valence sector it is a trivial 
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shift of the helds, 


Tyipv{x) = ipvix + jl), (A6) 


where /i denotes the unit vector in the /i direction. In the sea sector it involves 
transformations of the spin and taste degrees of freedom. Explicitly 


TP'Psiy) 

TfPshj) 


2 L 
1 

2 L 


-ff,5 ® ^5)'4’s{y) + (/ ® + 7/^5 ® ^5)'4’s{y + 2/i) , (A7) 

P’siy) + 7m5 ® C 5 ) + P’siy + 2h)(/ ® - 7/.5 8) 6) • 


Note that the translation operators for different directions commute when acting on 
staggered bilinears, so the order of them is irrelevant in the product in eq. (IA5jl . 

Applying the projection operator eq. (Hi results in many terms with derivatives, 
which we can neglect. Acting with Ty on the valence bilinear gives 

7pv{y)Tv'i/jv{y) ^ (^y(l/) + Vj?/>y( 2 /))rv'+ VjV’n( 2 /)), (AS) 


where we have introduced the usual nearest-neighbor forward difference operator in /i 
direction, Similarly, acting with Tg on sea quark bilinear produces many derivative 
terms. Using eqn. dSi and (Hi one straightforwardly establishes 


'^^si.y)'^s'4’s{y) ^ + derivative terms. (A9) 

y y 

Using these two results we hnd 


V 


04f(r5,rv^) 


04 f(f 5 , Ey) -|- derivative terms. 


(AlO) 


The matrix T 5 differs from T 5 only in the taste matrix: is replaced by the average 

1 

2 = 1 

where the sum runs over all sixteen elements of the Clifford algebra. Only the identity 
= I survives this average; for the other 15 taste matrices the average is zero. Thus 
we conclude that we only need to consider sea quark bilinears in eq. (IA4I1 that are 
taste singlets. 
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3. Next we impose the constraint that the operators mnst be color singlets. There are 
only two ways to form snch singlets. Either the color gronp generators in Tg and 
Ty are eqnal and a snmmation over the generator index is performed, or the identity 
matrix is inserted instead. 


4. Now we form linear combinations of (94f(r5,ry) that are singlets nnder the hyper- 
cubic symmetry group of the lattice (vr/2 rotations and reflections) and also charge 
conjugation. The transformation properties of the staggered helds in the hypercube 


notation are listed in Refs. 
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60| . Since the matrix T^ is trivial in taste space these 


transformations act in spin space only and their form is the same as for continuum 
Dirac spinors. We therefore hnd that the gamma matrices in T^ and Ty must be 
equal with their open indices being properly contracted in order to form scalars under 
rotations and reflections. 


5. Finally we select the operators that are invariant under the chiral symmetries. Each 
bilinear must be separately invariant under the full chiral symmetry group. This 
excludes all gamma matrices but the vector and axial vector. 


This procedure produces all mixed four-fermion operators without derivatives that are sin¬ 
glets under all lattice symmetries. Taking the continuum limit one hnds the four invariant 
operators listed in eq. ©• 
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